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Abstract 

We have computed the spectral number variances of an extended random 
matrix ensemble predicted by Guhr's supersymmetry formula, showing a non- 
monotone increase of the curves that arises from an " overshoot " of the two- 
level correlation function above unity. On the basis of the most general form of 
N- level joint distribution that meets sound probabilistic conditions on matrix 
spaces, the above characteristic may be attributed to the attractiveness of the 
pair potential in long range(£' > Thouless energy) of the underlying level gas. 
The approach of level dynamics indicates that the result is " anti-screening " 
of the level repulsion in short-range statistics of the usual random matrix 
prediction until the joint level distribution undergoes a phase transition (the 
Anderson transition). 

PACS numbers: 05.45.+b, 05.20.-y, 05.40.+j, 72.15.Rn, 02.50.-r 



1 



Typeset using REVTeX 



There have been considerable efforts in solid-state and random matrix theories(RMT) to 
formalize metal-insulater transition phenomena as regards the pertinent electron energy level 
statistics, which seek a powerful and unified method to generalize the standard Gaussian 
ensembles initiated by Wigner, Dyson and Mehta (see a comprehensive review on the 
recent development 0). 

Guhr has studied a method based on supersymmetry || , and recently obtained a formula 
for computing the two-level correlation function X 2 (r) for a Hamiltonian system H = Hq + 
aHi with a parameter A(A = A(a)) to describe a transition from regularity to chaos ||||; 
< A < oo. It is assumed that A = 0(also a = 0) represents the Poisson regularity for the 
system H alone (which is subject to the Poisson statistics), and A = oo the Gaussian unitary 
ensemble (GUE) for the system H\ alone (which is subject to the GUE statistics). Namely, 
in Refs. @,[|], a double integration representation of the X 2 (r, A) function is presented: the 
one involving a single Bessel function J\(z) || is shown to converge for A — > oo to the 
well-known GUS correlation function 1 — (siszei) anc i f or the lowest non-zero term of the 



X 2 (r, A) = — J exp (— k 2 /2J sm(rk/2)dk (valid for A ~ 0.1 or smaller). (1) 

Our concern in this report is not any inquiry about the supersymmetry basis of the 
formula nor its derivation, but a consequence of what Guhr called "overshoot" of the quantity 
X 2 (r, A) above unity, i.e. 

a range of r exists, where X 2 (r, A) > 1. (2) 

This can best be visualized by plotting X 2 (r, A) against r: here in Fig.l, we calculate four 
curves A = 0.1, 5, 10 and oo, and compare them with the corresponding curves which are 
provided by another formula of Gaudin's model due to us (Some feature of this model 
will be discussed later). 

Since the unfolded scale is used always in the present formulas and figures (which also 
ammounts to a change of the perturbation parameter a to A), inequality @ implies the 
negativeness of the corresponding cluster function Y 2 (r, A) |TJ so that 

Y 2 (r, A) = 1 — X 2 (r, A) < in the same range. (2') 

This causes a significant modification of the structure of the level statistics for such systems 
from the standard RMT. We would like to point out the following two remarks. 

1. Number variance curve based on the well-known formula 

S 2 ( S , A) = s - 2 f\ s - r)Y 2 (r, X)dr. (Y 2 (-r, A) = Y 2 (r, A).) (3) 

J 

It can be seen readily from Eq.fl3|) that ^-S 2 (s) = — Y 2 (s) < 0, which means that a S 2 (s) 
curve has an inflection point at the zero of Y 2 : more precisly, the increasing behavior of the 
curve changes its second derivative from minus to plus at the inflection point. 

2. N-level joint distribution which we assume to be of the form 



P(xi,x 2 , ..,x N ) = ^exp 



-(3 <t>{xj - x k ) 
j<k 



13 = 1, 2 and 4. (4) 
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Here, the sum J^jKk^^j ~ x k) in the right hand side of eq.(4) to represent the interaction 
between the levels is common in most of extended RMT nowadays, and we shall state its 
axiomatic basis in 1st aspect of Balian's strategy below. 

The cluster expansion theory for imperfect classical gases in statistical mechanics tells us 
that, if the binary potential for a pair of the molecules is positive (i.e. everywhere repul- 
sive), the second-order cluster function must be non- negative at least for the low density 
limit where it is approximated by the minus of the Mayer function, 1 - e _/3<?i(r) (/3 > 0) ( 
and a further account therein). Guhr's example of the smallest parameter value can 
be regarded as such a limiting situation so that the existence of the overshoot of Guhr's 
correlation function on the r-axis has a significance that, when the joint distribution of his 
model is expressed in the form of Eq.(f|), the potential <j) for a pair of levels Xj and Xk must 
have an attractive portion somewhere in the configuration space (xj,Xk) of the levels. Does 
this mean a level attraction rather than the level repulsion of the standard RMT sense? 

The possibility of long-range level attraction has been discussed by Jalabert, Pichard and 
Beenakker |8j for actual disordered metals. We shall discuss this question, referring to an- 
other recent paper by Weinmann and Pichard || who observed the behavior of non-monotone 
increase of the number-variance curves in a selected actual matrix ensemble and analyzed 
the data on their advocation of "Gaussian matrix ensemble with preferential basis" (cf. [TO]). 
Before going, we present a numerical and graphical confirmation of the non-monotone be- 
haviors of the theoretical £ 2 curves predicted by Guhr's formula (1) and that from a related 
approximation (see Eq. (4.3) in Ref. [PH ). 



As can be seen in Fig. 2 (a), the non- monotone characteristic is common for all the pa- 
rameter values, although the overshoot becomes obscure in Fig. 1(a) quickly as A increases. 
So, it is important to fix the range of A on which the non-monotone character remains to 
exist, in particular, to ask its existence for large A's. The latter question has been answered 



affirmatively by Fram, Guhr, and Miiller-Groeling |TTJ, who have provided a numerically 
tractable formula to replace Eq.(l) for X2(r, A). Moreover, they have shown a further sim- 
plification of the formula which turns out to be a divergent perturbation(i.e. starting from 
A = oo) by which £ 2 (s) for A ^> 1 can be understood easily. 

Leaving the actual graphs of confirmation to the paper ||11|| , we argue the existence of 



overshoot for A ^> 1 by means of their reduced formula: 

X 2 (r, A) = X° UE (r) + X^\r, A) + xj 2) (r, A) + (correction), 

where the term (correction) vanishes faster than r -2 for r — > oo. In terms of the cluster 
function Y"2(r, A) with neglect of such corrections, 

/ sin7rr\ 2 



y 2 (r,A)= ™ -X?\r,\)-X?\r,\), (5) 



where 



and 



\ 7IT / 



X?>(r,\) = — . 2 , 2j _ 2 , X^(r,A)dr = l, (5a) 

7r (7rA z r + r z ' 



1 r 2 - fvrA 2 ^ 2 
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It shows that the first-order correction X% (Breit-Wigner term) compensates the familiar 

(2) 

GUE part, together with the second-order correction X 2 contributing negatively to the 
cluster function for almost all r values(|r| > r with some small r ), which confirms the 
negativeness of Y 2 (r, A) above the r . At the same time, f*™ Y 2 (r, X)dr = for almost all A 
values, indicating that the E 2 curve approaches an asymptotic straight line that is the pure 
Poisson with no deviation of the coefficient from unity. The feature quite differs from the 
counter example of Gaudin's model |§, as exhibited in Fig. 2(b). 

To be a significance in the above two works, the authors of || and JIT] have presented 



actual matrix ensembles which display, as computer experimental results, the non-monotone 
number-variances: this stimulates us to construct an adequate statistical model capable of 
a unified description of both types; the monotone type and the non-monotone type. So, let 
us examine the concept of Gaussian matrix ensemble with preferential basis of Refs. |9||ll| • 
Our conclusion is precedently stated: the formulation as it stands is not capable of 
describing the non-monotone characteristic, because their potential function is everywhere 
positive, as exemplified by Eq. (|TT|) of Ref. [10| that is identical to || for Fig.2(b). However, 



this formulation can be revised without loss of sound mathematical basis to include the 
non-monotone characteristic, which we will discuss in the rest of the present report. 

There are two aspects of the Gaussian matrix ensemble with preferential basis to be 
reexamined seriously, in order to allow it with the possibility of attractiveness of the pair 
potential: the first is a proper introduction of the preferential basis, and the second, a right 
way of setting up constraints in maximum entropy principle for the Gaussian probability 
density function. Both aspects are combined to give the most efficient realization of Balian's 
strategy |T2] ]. 



1st aspect. An element of Balian's strategy in RMT is an introduction of Riemannian 
metrics on matix spaces in the form TidMdM* by which an information quantity is written 
(his postulate A). This is particularly relevant in the present problem, because a fixing of 
the metric tensor is the starting point of the whole subject of a Gaussian probability that 
is an exponential of the metric form. Namely, 

P(A) = Cexp[-l-K(A,A)}, K(A, A) = J2 K a,pA* a Ap > (equality onlyA = 0), (6) 

where A a is a tangent vector component of the metric form as the Gaussian random variable. 
We then ask what is the most general metric form defined on matrix spaces (See a detailed 
discussion [|13j). The answer is given by a covariant bilinear form of any two N x N 
hermitian(or, unitary) matrices A and B depending on another hermitian matrix H to 
satisfy representation invariance all together. It is expressed as 

K H (B,A)=ReTrB*C H (A), (7) 

in terms of a linear superoperator C# on A, and satiesfies the covariance condition 

Ku*hu(U*BU, U*AU) = K H (B, A) for any unitary U G Gp, (7a) 

where Gp is the symmetry group associated with f3. An adequate choice of the representation 
is naturally the if -diagonal representation ( H D = diag(xi, .., x^) by a choice of U to 
diagonalize H) that yields preferential basis: by this choice we can rewrite (7) as 
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K H (A,A) = J2 c ( x j)\ A jj\ 2 + 2 H f(xj, x k)\Ajk\ 2 with real positive c(x),f(x,y) (8) 

j j<k 

and f(x,y) = f(y,x). (8a) 

By definition of the trace operation in eq.(7), there is no actual priority of the repre- 
sentation basis in our "preferential basis", as it were supposed by the saying " if-diagonal 
representation" that might sound a kind of symmetry break[10] to prohibit some U G Gp. 
Instead, the unitary covariance (7a) allows U to cover the full group Gp so that the Gaussian 
variables {Aj^} are decorrelated from {xj}, and upon being integrated over these variables, 
the resulting reduced probability density function depends on the eigenvalue indices only 
through {xj} that are mutually fully equivallent to each other (so-called identically dis- 
tributed). Another characteristic of the Gaussian probability function (6), when the form 
(7) is inserted (with a = jk), is the statistical independence of A^s; the fact that the metric 
tensor in the form (8) is diagonal with respect to a. Let us further impose the third condi- 
tion of translational invariance(xj — > xj +x makes the metric tensor unchanged). Then, the 
form (8) is further simplified such that a single function f(x — y) is enough to characterize 
the iV-level joint probability function P{x\, ..,xn) as in the form (4). (This is seen by the 
integration of the function (6) over the cotangent variables defined by A^ = 2f(xj — Xk)Aj^.) 
We then get 

0(r) = --bg/(r) = -log-— y. (9) 

This is because the result of the integrations is expressed as the square-root of the determi- 
nant of the metric tensor that yields P(xi, .., xn). At the same time, a probabilistic meaning 
of the function f(xj — Xk) is assingned to be the variance of the Gaussian cotangent variable 
Ajk(oT, its each component, if it is complex, and quaternion depending on the multiplicity 
(3 = 2, and 4, respectively). To summarize, we can say 

the most general expression for N -joint level distributin Pn as an integral reduction of the 
Gaussian form (6) satisfying the identicalness, statistical independence and translational in- 
variance must be the form (4). (An inclusion of one-level potentials is another matter.) 
2nd aspect. It is important to remark that, in order to apply the maximum entropy prin- 
ciple to assign a form to the probability function P(A) (Balian's postulate B), a Gaussian 
probability has a special property of its entropy. Namely, the Gaussian probability func- 
tion Po(Ajk) with the variance f(xj — Xk) is characterized by the maximum entropy of all 
the probability functions with the same variance, thus denoting the entropy functional of 
P, (— \ogP)p, by H[P], we can express this fact as 

max H[P] = H[Pg] under constraint Var(Ajfc) = f(xj — Xk)- (10) 

(For the legitimacy of a maximum entropy principle with non-constant constraint, see [14].) 
It tells us that the variance (|H]) must be the right quantity of constraint for the present 
maximization problem. Therefore, without any further information about the matrix en- 
semble, we have no criterion about the degree of attractiveness of the pair potential <fi that 
is only related to the form of variance / as Eq.@: all what is needed is the positivity / > 0. 
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However, a powerful information can be provided by level dynamics ||(cf. [14]), where it 
is shown that, if the diagonalization process of a hermitian matrix H (or, its unitarization 
U^e ) is put into a Hamiltonian dynamics, a possible equilibrium state of this dynamical 
system can be selected by the above maximization as a one-parameter family of canonical 
distributions (in the sense of statistical mechanics of Hamiltonian systems) and represented 
by the form 



P G (A) oc [] exp 

j<k 

and hence 



1 - ,2 



k 



2 

jJiT 

with fir) = I -I (hermitian case), (11) 

1 + jjr 2 



llog|l + i. (Ha) 



2 fir 



A positive parameter value fx corresponds to Gaudin's linear-gas model [ 151 , also identical to 



the formula of Refs. PJT0|] for which is fully repulsive. To our emphasis, however, there is 
no a priori reason to restrict us to the positivity of fi (as far as f(r) is retained as positive), 
and let us look at the possibility of choosing a negative /i(This amounts to a sign change of 
Yukawa's parameter 7 |16[) : it makes the potential partially attractive such that 



11 / — 11 / — 
0(r) = -l g(— — - 1) \r\<l/J\n\', — log(l — -j— j— 2) \ r \>l/yM> 
2 \jj,\r z v 2 |/i|r z v 

with the attractive range 1/^2|//| < \r\ < 00, (12) 

and at the same time, for small values of the level distance including the repulsive range(|r| < 

l/\/2H), 

/(r)=r 2 |/i|/(l-|^k 2 ) \r\<l/y^\. (13) 

(A schematic behavior of the entire potential function 0(r) can be seen in Fig. 3.) This 
can be interpreted as an anti- screening of the level repulsion in short-range rather than the 
screening assigned to it in Gaudin's model by the original proposal flOf , and as the actual 



level attraction in the range specified by (12) for long-range level statistics. 

There exist two routes of the Wigner-Dyson Gaussian standard of electron energy level 
statistics to Poisson regularity; one /i(= 0) tends to infinity on the positive axis, and the other 
/i(= 0) to infinity on the negative axis, and it is remarkable that the latter route undergoes 
a singularity of the level distribution that can well be assigned to be the Anderson transition 
with Thouless energy E c = A/ \f\/~i\(A is the average spacing in energy scale ): the resulting 
Poisson statistics must be on localized electrons that, at the same time, give rise to a real 
symmetry-breaking of the representation (the metric tensor then depends on the indices jk 
as e~ a '- ?_fc '). Existence of both repulsive and attractive inter-level pair- wise interaction was 
first indicated by || for disordered metals, and from the present viewpoint these provide 
a distinction between absence and presence of the transition, and it can be tested by the 
monotonicity or non-monotonicity in the number-variance curve as in Fig.2. 
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FIGURES 



Fig 1 
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FIG. 1. X2(r, A) for different values of the transition parameter A from two different theories, 
(a) Calculated from Guhr's formular for A 1 and ir\ 2 » 1. (b) Calculated from Gaudin's model 
(Eqs. (4.5) of Ref. [6] ) From the left to right the curves corresponde to A = 0.1, 5, 10 and oo. The 
curves of A = 5 and 10 are almost indistinguishable from the GUE curve (solid line). 
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Fig 2 
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FIG. 2. Number variance S 2 (s, A) for two different models, (a) From Guhr's model (Eq. (4.4) 
of Ref. [10]). (b) From Gaudin's model due to Hasegawa and Ma [6]. The non-monotonic behavior 
in (a) is very obvious. From the left to right the curves corresponde to the Poisson, A = 0.1, 5, 10 
and oo (GUE), respectively. 



9 




FIG. 3. <fi(r) vs r for ji = —0.5 given by Eq. (11a). 
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